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ABSTRACT 

The Australia Telescope Compact Array (ATCA) has been used to search for methanol 
maser emission at 6.7 GHz towards OMC-1. Two features peaking at 7.2 km s _1 and 
— 1.1 km s _1 have been detected. The former has at least two components close in both ve- 
locity and position. It is located south-east of the Orion Kleinmann-Low (Orion-KL) nebula 
in the region of outflow traced by the 25-GHz methanol masers and the 95 -GHz methanol 
emission. It is shown by modelling that in contrast to the widespread opinion that simulta- 
neous masing of methanol transitions of different classes is impossible there are conditions 
for which simultaneous masing of the class II transition at 6.7-GHz and some class I tran- 
sitions (e.g. the series at 25 GHz) is possible. A relevant example is provided, in which the 
pumping occurs via the first torsionally excited state and is driven by radiation of the dust 
intermixed with the gas in the cloud. In this regime the dust temperature is significantly lower 
(T«60 K) than in the case of bright 6.7-GHz masers (T>150 K). The narrow spectral feature 
at —1.1 km s _1 has a brightness temperature greater than about 1400 K, which suggests that 
it is probably a maser. It emanates from the Orion South region and is probably associated 
with the approaching part of outflow seen in CO. The 25-GHz maser associated with OMC- 
1 was observed quasi-simultaneously with the 6.7-GHz observations. No 25-GHz emission 
associated with the —1.1 km s _1 6.7 GHz feature towards Orion South was detected. 
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1 INTRODUCTION 

The Orion Molecular Cloud 1 (OMC-1) is a well known region 
of massive star formation containing several young stellar objects 
(YSOs) at the very early stages of their formati on as well as out- 
flows related to them (see, e.g. the review by iGenzel & Stutzkil 
1989). Methanol masers in a variety of transitions are often found 
in such regions. Masing transitions are traditionally divided into 
two classes. Although masers of both classes often co-exist in the 
same star-forming region, they are usually seen apart from each 
other I Men ten 1991aJ3). 

OMC-1 harbours strong class I methanol masers. These 
masers are believed to trace distant parts of the outflows from YSOs 
and are ass ociated with the early stages of the massive star for- 
mation (e. g. . iMentenl 1 99 hk lEllin gsen 2005). The strongest (about 
150 Jy) masers of the J2— Ji E line series near 25 GHz are ob- 
served in OMC-1 and have inspired numer ous ob servational and 
theoretical studies (see, e. g., |MenterUgUil,ir i988ct Ijohnston et alj 
ll992HSobolev. Watson & Okorokovl2003l) . 



E-mail: Maxim.Voronkov@csiro.au 



In contrast, strong class II masers have not been detected 
in OMC-1. Only weak emission in the strongest class II 5i — 
60 A + transition at 6.7 GHz has been reported from this source 
bv lCaswell et alj ll995t) . and it was unclear whether the observed 
spectral line is a maser or has a quasi-thermal origin (i.e., the level 
populations are without inversion). Currently, class II methanol 
masers are thought to appear in the vicinity of the massive young 
stellar objects (YSOs) at the e arly stages of ultra-com pact HII re- 
gion tIJCIirrj_jyohition (e.g., lEllingsen et alJll99(J : iMinier et all 
l200ltlCodella et all2004l) . The transition at 6.7 GHz manifests flux 
densities which are the highest among methanol masers and can 
reach 5000 Jy. Maser emission from this transition is widespread 
and to date has been detected tow ards more than 5 00 sources 
(see, e.g. the compilation bv lMalvshev & Sobolevl2003i) . Being the 
strongest this transition is the best prospect for searches for new 
class II methanol maser sites. 

Evolutionary stages of the massive YSOs at which differ- 
ent classes of m ethanol masers exist do intersect considerably. 
lEllin gsen (2005) detected class I methanol masers at 95.1 GHz to- 
wards 38% of a statistically complete sample of the class II masers 
at 6.7 GHz. No extensive search of class II masers has been per- 
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Table 1. Dates of observations and a summary of array configurations. The 
compact H75B configuration featured the north-south spur of the ATCA 
providing a better uv-coverage for equatorial sources. In the 1.5D configu- 
ration session the CA05 antenna did not observe and, therefore, the second 
shortest baseline CA04-CA05 was not available. Other sessions had all 6 
antennae working. 



Date 


Array 




Baseline length (m) 


UT 


configuration 


Min 


Max 


Max 








no CA06 


with CA06 


9- 10 June 2003 


H75B 


31 


89 


4408 


12 November 2003 


1.5D 


107 


1439 


4439 


5 January 2004 


6A 


337 


2923 


5939 


10 April 2004 


EW367 


46 


367 


4408 



formed towards the positions of isolated class I masers. This is 
partly due to the lack of information about such sources since most 
class I masers have been found towards known clas s II masers, 
which are often located within a single dish beam (e.g.. lSlvsh et alJ 
11994 . 

Therefore, among known sources OMC-1 is somewhat ex- 
ceptio nal: it posesses stron g class I masers produced by an out- 
flow lJohnstonetai] Fl992). while there are no strong class II 
masers associated with the YSO producing this outflow. Re- 
cent observations have shown that some hydroxyl masers and 
class II methanol masers can als o be associated wit h out- 
flows lArson. Reid & MenterJ 120031: ICodella et all [2004). This 
is supported by statistical analysis of correlation between the 
mas er velocities and those o f the molecular outflow trac- 
ers (Malvshev & Sobolev 2003) and is reasonable because the 
methanol abundance is significantly increased in the interface re- 
gions, where the outflow interacts with the amb ient material (e.g., 
Bachi ller et all 1 99.41 1 99?3:IOibb & DaviJl998h . It is presently un- 
clear whether weak class II masers can co-exist with the class I 
masers or they always trace different parts of outflow. We have un- 
dertaken a sensitive search at 6.7 GHz in order to figure out whether 
there are weak class II masers in the vicinity of the YSOs and out- 
flows in OMC-1. The 25-GHz class I maser in OMC-1 was reob- 
served to check its current status. 



2 OBSERVATIONS 

The observations were made using the Australia Telescope Com- 
pact Array (ATCA). Details on the array configurations 1 and ob- 
serving dates are summarized in Table Q During the first ses- 
sion on 2003 June 9-10, the 25-GHz maser was observed quasi- 
simultaneously with the 6.7-GHz observations in 7 short scans 
switching the frequency approximately every 90 minutes. The pri- 
mary aim of the subsequent sessions was to estimate the 6.7- 
GHz source size and therefore the 25-GHz transition was not 
observed. The adopted rest frequencies are 24959.080 MHz and 
6668.5192 MHz fo r 5 2 - 5i E and 5i - 6o A+ methanol tran - 
sitions respectively lLovasll992l:lBreckenridge & Kukolicrl 19951) . 
For both frequencies and all sessions the correlator was configured 
to split the 4-MHz bandwidth into 1024 spectral channels yield- 
ing a spectral resolution of 3.9 kHz or 0.05 km s _1 at 25 GHz 



and 0.18 km s _1 at 6.7 GHz (no Hanning smoothing has been ap- 
plied). Two orthogonal linear polarizations were observed at each 
frequency, except the CA01 antenna in the first session on 2003 
June 9-10 which observed only horizontal polarization (at both 
frequencies) due to technical problems. During the data reduction 
the source has been assumed to be unpolarized and both polar- 
izations, if presented, were averaged together. The position of the 
phase and pointing centre was 02000 = 5 h 35 m 15=069, £2000 = 
-5°23'42'.'53 for the first session and 02000 = 5 h 35 m 15 s 439, 
^2000 = -5°23'07'/0 for the latter observations. The full width at 
half maximum (FWHM) of the primary beam was 7!2 at 6.7 GHz 
and 2!0 at 25 GHz. The global pointing solution was used for both 
frequencies giving an rms pointing accuracy of 5 — 10". 

Most of the data reduction was performed using the MIRIAD 
package (13 May 2004 release) following standard procedures with 
the exception of bandpass calibration for the 6.7-GHz observations. 
We found that the accuracy of the bandpass calibration limited our 
ability to detect a weak spectral line. The shape of the uncali- 
brated bandpass was reasonably straight within the small part of 
the band occupied by the source spectrum and it was therefore safe 
to skip the bandpass calibration. However, this approach fails on 
the shortest baselines. The continuum emission is very strong to- 
wards OMC-1 at short spacings. Therefore, it is difficult to subtract 
this continuum accurately without a proper bandpass calibration, 
which would require unrealistically long observations of the band- 
pass calibrator. A residual error in the continuum subtraction tends 
to create weak and broad spurious lines in the spectrum. A similar 
problem makes it extremely difficult to observe a weak spectral line 
towards a strong continuum source using a single dish. Therefore, 
all the 6.7-GHz data resulting from baselines with a projected spac- 
ing shorter than 50m were flagged. This has a serious impact on the 
6.7-GHz observations using the H75B array configuration, where 
about one half of the data was discarded. Shadowing was also a 
problem in the compact H75B array configuration throughout the 
whole session. All the 25-GHz data from baselines with a shad- 
owed antenna have been flagged (3 per cent of visibilities from the 
total number of baselines without the CA06 antenna). The rejec- 
tion of short spacings at 6.7 GHz discards all shadowed baselines 
automatically. During imaging, natural weighting was used to min- 
imize the noise in the image and to provide better stability in the de- 
convolution algorithm (CLEAN) for the near-equatorial source. The 
data resulting from the long baselines with the CA06 antenna were 
used for non-imaging analysis only. In all sessions except the first 
one the standard ATCA primary calibrator 1934-638 was used for 
the amplitude calibration. For the first session 0605-085 with the 
adopted flux density of 1.8 Jy at 6.7 GHz was used instead due to 
the siderial time range allocated. The 25-GHz data were corrected 
for the atmospheric opacity using archival meteorological data and 
the flux density scale was calibrated using observations of Uranus. 
The PLBOOT task of MIRIAD was used to compute the scale fac- 
tor from the planetary observations using a built-in ephemeris. For 
the date of 25-GHz observations it implied the total flux density of 
Uranus equal to 0.76 Jy. The continuum source 0539-057 was used 
as a secondary calibrator at both frequencies. Although the calibra- 
tion and imaging were performed in the MIRIAD, the AIPS++ pack- 
age developed by NRAO 2 was used for a direct manipulation of the 



1 More details on the ATCA configurations and base- 
line length distribution are available on the web 

http://wwwMarrabriMtnfxsiroMu/opemtions/army.configurations/configuratu 



2 The National Radio Astronomy Observatory (NRAO) is a facility of the 
National Science Foundation operated under cooperative agreement by As- 
sociated Universities, Inc. 
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Figure 1. A cumulative image of the 6.7-GHz methanol emission in OMC-1 (all data except those resulted from the CA06 antenna have been used) and 
the spectra (a vector averaging of all data corresponding to uv-distances less than 5 kA has been used) at the peak positions for each feature. The filled 
ellipse represent the FWHM of the synthesized beam of 27'.'4x lff'2 and its position angle of 0?8. Contours are 20, 30, 40, 50, 60, 70, 80 and 90 per cent of 
0.49 Jy beam -1 . The two filled circles represent IRc2 (lower) and BN (upper) objects. 



Table 2. Results of the Gaussian fitting to 6.7-GHz images of individual spectral channels. All uncertainties are formal errors of the Gaussian fit. A good 
estimate of the actual lcr positional accuracy is 5" x2" (full width) with a positional angle of 0?8 (the synthesized beam has the FWHM of 27.4" X 10.2" and 
the same positional angle). 
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22'49"±6 
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22'42"±2 


0.22±0.02 


0.57 


43±8 


17±3 


0±6 


33±10 


14±4 


— 1=1=21 


7.9 


14.37i0.06 


22'46"±2 


0.17±0.01 


0.42 


40±8 


17±4 - 


5±8 


29±11 


13=1=5 


-10±32 



6.7-GHz data to investigate the behavior of the visibility amplitude 
at different uv-distances. 

3 RESULTS 

3.1 The 6.7-GHz data 

To analyze the spatial and velocity distribution of the 6.7-GHz 
emission in OMC-1 we constructed a spectral cube and vector av- 
eraged spectra at different positions. All data from all sessions, ex- 
cept those including the CA06 antenna were used to produce the 
cube. A vector averaging (i.e. averaging of complex visibilities with 
phase) was used for the spectra because the signal was not strong 
enough for other methods of presentation. To avoid decorrelation in 



the spectra only data corresponding to uv-distances less than 5 kA 
were included. A cumulative image of the source obtained from 
the cube is shown in Fig. Q along with the spectra. The brightness 
of each pixel is the maximum across the velocity axis (different 
planes of the cube correspond to different spectral channels and, 
thus, to different velocities). The image in Fig.Qshows that there 
are two sources separated by about 80" in the north-south direc- 
tion. Given the positions of the sources, the reduction in amplitude 
due to primary beam attenuation is estimated to be less than 7 per 
cent, which is less than the precision of the absolute flux density 
calibration. This effect can therefore be neglected at 6.7 GHz and 
no primary beam correction has been applied to either the spec- 
tra or the image. These sources in the map in Fig. Q correspond 
to two spectral features, a broad one peaking at 7.2 km s" 1 and 
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a narrow one peaking at —1.1 km s _1 . An inspection of individ- 
ual planes of the spectral cube confirms that no emission brighter 
than 0. 15 Jy beam -1 is visible at the position of one feature within 
the velocity range of the other. Both features were detected in the 
vector averaged spectra constructed separately for each individual 
observing session, which can be considered as an independent con- 
firmation. Radial velocities in the range 4—9 km s _1 are typical 
for the Orion region. The emission from the hot core component 
peaks in the range 3—8 km s -1 , while the me thanol emission core 
(compact r idge) pea ks around 7—9 km s~ - 1 lMentenetal]|l986t 
Wrigh t et alll996tlCragg et all200ll) . lMenten et alHl988bl) found 
the methanol emission in the southern source peaking at approxi- 
mately 6 km s . However, no methanol emission has been previ- 
ously reported at or near — 1 . 1 km s ~ 1 . 

The northern source (the feature peaking at 7.2 km s _1 ) is lo- 
cated about 10" to the south of the IRc2 object. This position is 
typical for methanol l ines at these velocities and is known as the 
compact ridge I Ment en et"allll988rlcl IWrigfit et alJll996h . From 
Fig Qthe source appears extended and two components are readily 
distinguished. The spectrum also supports this as the profile of the 
broad feature is asymmetric. To reveal the velocity structure of the 
source we fitted a Gaussian source model to individual planes of 
the cube used to construct a cumulative map shown in Fig.Q The 
results of these fits are listed in Table|2| The first five columns are 
the LSR velocity of the channel, the position and the peak flux den- 
sity of the fitted Gaussian and the total flux density of the source. 
The last six columns are the Gaussian parameters, the length of the 
major and minor axes (at half maximum) and the position angle of 
the major axis, for the fitted Gaussian and the result of deconvolu- 
tion of the synthesized beam. The uncertainties listed in the table 
are the formal errors from the Gaussuan fit. The formal errors for 
the deconvolved Gaussians are obtained under assumption that the 
synthesized beam is known precisely. A good estimate of the actual 
lcr positional accuracy is the synthesized beam size divided by the 
signal to noise ratio in the image (e.g.. lFomalonT l999). The noise 
level in the map of each spectral channel (cumulative map is shown 
in Fig.0 is about 20 per cent of maximum. The FWHM of the syn- 
thesized beam is 27" 4 x \$!2 and the position angle is 0?8, which 
gives an estimate of the lcr positional accuracy of 5" x 2" with the 
same position angle. Note, however, that for a blended component 
like 7—8 km s" 1 feature is, the actual accuracy may be worse. The 
fitted positions in the Table|2|show a trend of moving towards the 
east as the velocity decreases from 8 km s _1 to 7 km s _1 . This 
supports the idea that the source is a blend of two components with 
similar velocities and positions. The southern source (—1.1 km s _1 
feature, Fig.0 is located at the position of the Orion-South region. 
This re gion is a known source of methanol emission l Ment en et alJ 
1988b), however, the thermal emission peaks at a velocity of ap- 
proximately 6 km s _1 . The results of the Gaussian fit into the im- 
ages of different spectral channels across the line profile of the 
— 1.1 km s" 1 feature are the same within the fitting uncertainty. 
Therefore, in the Table|2|we listed only one of the fits, correspond- 
ing to the peak velocity. 

A direct analysis of the visibility data is another, complemen- 
tary way to obtain information about the source structure. Given a 
model of the source, this method is able to produce more reliable 
uncertainties than fitting a deconvolved image. For a weak source, 
however, one is limited to use the simplest source models only as 
the ambiguity of the fit increases rapidly with the number of free 
model parameters. Fig.|5|contains the dependencies of the visibil- 
ity amplitude on the baseline length for the 7.2 km s _1 (upper plot) 
and —1.1 km s _1 features (lower plot). Each point in the plot rep- 
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» EW367 




2x10" 4x10" 6x10" 8xl0 4 

Baseline length, A 



Figure 2. The dependence of the visibility amplitude on the baseline length 
for 7.2 km s — 1 (upper plot) and —1.1 km s — 1 (lower plot) 6.7-GHz fea- 
tures. Different symbols correspond to the data from different array config- 
urations. The solid and dashed lines represent the results of model fitting. 
The EW367 data for the —1.1 km s - 1 feature were fitted separately (dashed 
lines) from the data from other configurations (solid line). 



resents a single baseline. Baselines from the different array config- 
urations are shown by different symbols. All visibility data were 
vector averaged together separately for each baseline using the po- 
sitions shown in Fig.Qfor the 7.2 km s _1 and —1.1 km s -1 features 
respectively, regardless of the projected length. The horizontal er- 
ror bars in Fig. [5] represent the scatter of the projected length for 
each baseline. The vertical error bars show the scatter of the aver- 
aged amplitude. This complicated scheme of averaging is dictated 
by the weakness of the source. Different fit results are shown by 
the solid and dashed lines in Fig.|2| The 7.2 km s -1 feature can be 
described by a symmetric Gaussian source model with a total flux 
density of 0.6±0.1 Jy and FWHM of 18" ±9". Taking into con- 
sideration that the size of this symmetric model is an average size 
as the actual source is most probably asymmetric, the value is in 
agreement with Table|2|(deconvolved sizes and separation between 
peaks at different velocities). From Fig.|2|the —1.1 km s _1 feature 
appears to have undergone a flare during the EW367 observations 
(2004 April 10) as the flux densities are significantly different from 
other configurations (including the 6A observations from 2004 Jan- 
uary 5). The solid line represents a point source model with a flux 
density of 0.5±0.1 Jy fitted to all data except that from the EW367 
configuration. One can see in Fig.[2]that the —1.1 km s _1 feature 
is unresolved for baseline lengths of about 8 x 10 4 A, which im- 
plies the source size less than about 2" 6. The EW367 data were 
fitted separately. A gradual decrease of the visibility amplitude at 
the short EW367 baselines seen in Fig.|2|encouraged us to try a two 
component model, although the point source model is adequate for 
the data (with the observed noise level an extended component is 
not necessary to explain the plot). The point source model (short 
dashed line in Fig. |5J implies the flux density of 1.2±0.2 Jy. The 
two component model (long dashed line) includes an unresolved 
component with the flux density of 1.0±0.2 Jy and a symmetric 
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Figure 3. Images of the 25-GHz maser (dashed contours) in OMC-1 at different velocities overlayed on the images of the 95-GHz emission (solid contours) 
obtained by Wrig ht et all ll99(Ji . courtesy of ADIL. The filled ellipses in the bottom right corner represent the synthesized beams at 25 GHz (large ellipse, 
28'.'7x 18'.'1, -64°) and 95 GHz (small ellipse, 7.1x4.9, 36°). Contours are 10, 20, 30, 40, 50, 60, 70, 80 and 90 per cent of 86.8 Jy beam" 1 at 25 GHz 
and 20, 30, 40, 50, 60, 70 80 an d 90 per cent of 40.4 Jy beam -1 at 95 GHz. The spectral resolution is 0.5 km s -1 and 0.05 km s _1 for the 95-GHz data 
obtained by Wright et "aTl ll996}l and our 25-GHz data, respectively. Bold open ellipses represent the 3<r uncertainty of the 6.7-GHz emission position and 
co rrespond to the fit resu lts at 7.2 km s _1 and 7.5 km s _1 for the left and right ellipse respectively. The filled circles represent 25-GHz maser spots revealed 
by Johnston et al. 1 19921) . For each subplot, only spots with the peak velocity closer to the subplot velocity than their half width at half maximum plus the 
width of one VLA spectral channel (0.15 km s _1 ) are displayed. The circles of a smaller radius correspond to faint 25-GHz spots (less than 10 Jy of the total 
flux). 



Table 3. Results of Gaussian fits to 25-GHz images of characteristic channels shown in Fig|3] All uncertainties are formal errors of the Gaussian fit. A good 
estimate of the actual \a positional accuracy is 2.4" X 1.5" (full width) with a positional angle of —64° (the synthesized beam has the FWHM of 28.7" X 18.1" 
and the same positional angle). 
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7.34 


14.87±0.03 


44.0±0.4 


35±2 


37 


29.6±1.1 


18.6±0.7 
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4±3 
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14.83±0.03 
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51 


30.4±1.2 


18.4±0.7 
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10±4 
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7.95 
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— 51=1=14 


8.28 


14.02±0.04 
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89±4 
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-53d=2 


28±2 


4±4 
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31.8±0.3 
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35.6±0.9 
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22±1 
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-35±7 


9.92 
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14=1=1 


17 


29.9±1.3 
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7±5 


-5±67 



Gaussian source with a total flux density of 0.4±0.1 Jy and the 
FWHM 17" ±13". Note, that the flux densities and source sizes in 
these models are overestimated (less for brighter signals) due to 
amplitude bias. However, the magnitude of this systematic offset is 
well below the random noise. 



3.2 The 25-GHz data 

The 25-GHz methanol maser in OMC-1 w as the first maser 
from the methanol molecule de t ected in space iBarrett et alll97ll 
IChui et alll9 7l. lMenten eTaTl fl988c) observed several lines (J= 
2 to 8) in the series and found that maser emission dominates in 
the spectra for the J= 4 line and higher while most of the emis- 
sion in the J= 2 and 3 transitions is of a thermal (non-maser) na- 
ture. Emission was detected from both Orion-KL and Orion-South, 
resembling the morphology of the source seen in mm-wavelength 
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methanol lines jMenten et allll988bl) . Ijohnston et alJ < 19921) im- 
aged the J= 5 and 6 masers in the KL region using the VLA D- 
configuration and showed that the maser consisted of a large num- 
ber of spots with a small velocity spread, distributed in a 40" long 
crescent shaped region. The spatial resolution of the ATCA in the 
H75B configuration is much lower than that of the VLA and is 
comparable with the resolu tion of the single dish Effelsberg map. 
Therefore, unlike Ijohnston et"al]l 19921) we cannot disentangle any 
information about individual maser spots and here only attempt to 
analyze a general picture. From the shape of the spectra on differ- 
ent baselines we have chosen several spectral channels correspond- 
ing to the peaks of distinct features. Images of the spectral chan- 
nels at velocities of 7.34, 7.58, 7.95, 8.28, 8.75, and 9.92 km s" 1 
are shown in Fig. [3] All these images have been co rrected for pri- 
mary b eam attenuation. The maser spots revealed bv ljohnston et alJ 
119921) at close velocities are overlayed on our 25-GHz image in 
Fig.|3| For each subplot, we show only those spots, which have peak 
velocities closer to the s ubplot velocity than th e half width at half 
maximum determined by Johnston et al. < 19921) for these spo ts plus 
the width of one spectral channel in the Johnston et al. 1 1992) data 
(0.15 km s _1 ). InFig.|3we also show the 95-GHz (the 8 - 7i A+ 
methanol transition) images sam pled at subplo t velocities from the 
spectral cube obtained by Wrig ht et alJ 1 19961) . which is available 
from the NCSA Astronomy Digital Image Library (ADIL) 3 and has 
a spectral resolution of 0.5 km s~ 1 (our 25-GHz data have a spec- 
tral resolution of 0.05 km s _1 ). In a similar manner to the 6.7-GHz 
data we fitted a Gaussian source model into each spectral channel. 
The results are listed in Table [3] The columns are the same as for 
Table [2] As before the uncertainties in the table are formal errors 
of the Gaussian fit and a good estimate of the actual la accuracy is 
given by the synthesized beam size and signal to noise ratio, which 
gives 2.4" x 1.5" with a positional angle of —64°. 

Although we know from the image of l.Tohnstonetal.Hl992l) 
that the 25-GHz emission contains a larger number of components, 
these characteristic velocities we show are sufficient to analyze the 
morphology o f the source within th e spatial resolution provided 
by the ATCA. iMenten et alJ fi988ct) also found 6 components in 
their Effelsberg data, which has a comparable spatial resolution 
to our measurement. However, the component velocities listed by 
iMentenetaj ] 11988 c) are slightly different. Moreover, our data im- 
ply the presence of a 7.58 km s _1 feature, which is not listed by 
Ment enet alJ ll988ch . but do not require a 9.6 km s _1 feature de- 
tected by them. The difference is most probably due to the interfer- 
ometric nature of our measurement as the output signal for a com- 
plex source depends on the size of individual baseline and its ori- 
entation with respect to the source. Another plausible explanation 
of the difference is te mporal variability in the source. The image of 
Johns ton et al] IY992) does not contain a spot at or near 9.6 km s _1 . 

A comparison of o ur 25-GHz image with the maser spots 
of Ijohnston et al J < 19921) in Fig. |5| (see also Table [5} shows that 
they are in general agreement. However, the south-eastern group 
of maser spots dominates in our image at low velocities. The most 
north-western spot in the 7.34-km s _1 subplot in Fig.|3|is, there- 
fo re, w eak according our data. This is the spot number 1 1 in Table 2 
of Ijohnston et alJ ll992h . This spot is the brightest spot observed 
(total flux 92 Jy) and is rather broad in the velocity domain (peak 
is at 8.0 km s _1 , half peak linewidth is 1.03 km s _1 ). Therefore, 
we suspect that significant temporal varia tions have occurred in the 
25-GHz masers since the observations of Ijohnston et alJ jl992t) . 
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Figure 4. The spectrum scaled to correct for the primary beam attenuation 
and the spectrum of the image flux density (bold line) of the 25-GHz maser 
in OMC-1. The subplot shows the 25-GHz spectrum at the position of the 
6.7-GHz emission near —1 km s — 1 . The scaling factors are 2.71 for the 
main plot and 1.16 for the subplot. 



The vector averaged spectrum of the 25-GHz emission is 
shown in Fig. The data on all baselines except those including 
the CA06 antenna, were vector averaged together for the posi- 
tion 02000 = 5 h 35 m 14^02, Jaooo = -5°22'3079, which corre- 
sponds to the peak of the 25-GHz emission in the image (Fig. [5} 
at 8.28 km s _1 . The amplitude was scaled by a factor of 2.71 to 
account for primary beam attenuation at the offset of the above 
mentioned position with respect to the pointing centre. However, 
this method of correcting for the primary beam is valid only for 
a point source. In a complex source the spots, which lie closer to 
the pointing centre than the above mentioned position, will appear 
brighter, while the effect of the spots which lie further will be un- 
derestimated. Moreover, the measured flux density is very sensitive 
to the pointing errors when the source is offset from the pointing 
centre. The magnitude of this effect is determined by the derivative 
of the primary beam and results in an error of 17—39 per cent (at 
the position of the 25-GHz peak) if a pointing error of 5 — 10" is 
assumed. The impact of the pointing errors on the measured flux 
density for a complex extended source can be alleviated by proper 
mosaicing observations. Another way to obtain a spectrum of the 
source from interferometric data is to sum the flux density in the de- 
convolved image for each spectral channel. This enables a proper 
primary beam correction to be done, taking into account the differ- 
ent positions of individual features and produces a single dish-like 
spectrum where the individual spectral features are emphasized. 
However, this method includes a non-linear deconvolution algo- 
rithm, which makes it difficult to analyze the noise and to assess 
whether a feature in the spectrum is reliable. In addition, in the 
case of a poor uv-coverage the quality of such spectrum is limited 
by the dynamic range of the image (5:1 in our 25-GHz measure- 
ment) and the absolute noise level is different at the line centre and 
at the wings. Freedom in the choice of the image region where pixel 
values are summed and the fact that the values of adjacent pixels 
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usually correlate make this method inappropriate for weak signals 
like the 6.7-GHz line discussed in the previous section. However, 
such a spectrum can be constructed for the 25-GHz maser which 
is bright. The spectrum obtained by summing over a 100" square 
centered at q 20 oo = 5 h 35 m 14=33, <5 20 oo = -5°23'14'.'6 is shown 
in Fig.|4]by the bold line. The size and position of the image region 
involved in summation is chosen to accomodate all the emission 
seen in Fig|5|and to avoid as far as possible the region near the edge 
of the field of view where the noise of deconvolution algorithm is 
amplified by the primary beam correction procedure. The most sig- 
nificant difference between this spectrum and the vector averaged 
one is near 10 km s _1 . This is because the source is offset to the 
north at these velocities (see Fig. [5} and the primary beam correc- 
tion uses the position of the peak at 8.28 k m s~ (see the discus- 
sion above). The spectrum resembles that of Johnston et al. 1 1992) 
obtained using the same me thod, although our sp ectrum peaks at 
8.28 km s" 1 , while that of Ijohnston et alj 1 19921) peaks at about 
7.9 km s _1 , and an additional feature near 7.6 km s _1 seems to be 
present in our spectrum. These differences can be explained by both 
the errors resulting from the dynamic range and the pointing inac- 
curacy as well as by different uv-coverage in the two experiments. 
Temporal variability of the maser is another possible explanation. 
The peak flux density determined using both the vector averaged 
spectrum and the image fl ux spe ctru m is in agreement with obser- 
vations of Johns ton et alj Il992r) and lMenten et alj (1988c) within 
the accuracy of our measu rement. 

iMenten et alJll98 8c1 report 25-GHz methanol emission from 
the Orion-South region (southern 6.7-GHz source) for the J= 2, 3 
and 4 lines in the series. The source has also mm-wave length 
methanol line emission peaking at 6 km s _1 (Menten et al. 1988c). 
However, no methanol emission at —1.1 km s _1 has been reported. 
The subplot in Fig. |4| shows the 25-GHz vector averaged spec- 
trum at velocities near —1.1 km s _1 . As for the main plot, all data 
except those resulting from the baselines with the CA06 antenna 
were vector averaged together. However, for the subplot in Fig.|4| 
the position of the averaging centre was 0:2000 = 5 h 35 m 13=93, 
$2000 = —5°24'05'.'3, which correspond to the position determined 
for the 6.7-GHz emission at —1.1 km s _1 (Fig. Q. The pointing 
centre was much closer to this position and, therefore, a smaller 
primary beam correction factor of 1.16 has been applied. The flux 
density of the 52 — 5i E emission at 25 GHz near the velocities of 
— 1.1 km s _1 is below the detection limit. The rms noise (la) is 
0.12 Jy. 



4 DISCUSSION 

4.1 Temporal variability and previous observations 

The non-imaging analysis discussed in 13.11 shows that the 
— 1.1 km s~ feature was stronger in the EW367 data (2004 April 
10) with respect to the data from other array configurations. This 
difference is unlikely to be due to instrumental effects. The derived 
secondary calibrator flux densities at 6.7-GHz are 0.68 Jy, 0.87 Jy, 
0.87 Jy, and 1 .07 Jy for the 2003 June 9-10 (H75B), 2003 Novem- 
ber 12 (1.5D), 2004 January 5 (6A), and 2004 April 10 (EW367) 
observations, respectively. These numbers imply that the precision 
of the absolute calibration is better than 16 per cent. Although the 
secondary calibrator appears brighter in the EW367 data, the differ- 
ence in the flux density is not enough to explain the —1.1 km s _1 
feature being 2—3 times stronger in these data. In addition, there is 
a good agreement between the amplitudes on short and long base- 
lines of the same session (Fig.|2j- Hence, a difference in the phase 



stability at different dates is unlikely to explain the observed dif- 
ference of amplitudes (a vector averaging of the data is used before 
the amplitudes are calculated) because the decorrelation should be 
stronger on longer baselines. A final argument against instrumen- 
tal effects is an agreement between the 7.2 km s _1 feature data 
from the different array configurations. Therefore, it is likely that 
the —1.1 km s" 1 feature had a flare in April 2004. 

The fi rst observation s of th e OMC-1 at 6.7 GHz were car- 
ried out bv ICaswell et alJ 1 19951) on 1992 December 17-24 us- 
ing the Parkes 64-m telescope. The telescope was pointed towards 
0:2000 = 5 h 35 m 14 s 35, S 2a00 = -5°22'24'/4. A weak (0.5 Jy) 
feature at 7 km s~ was re ported. The detection limit was about 
0.3 Jy <Caswelletallll995h . So the flux density measured more 
than a decade ago at Parkes for the 7.2 km s _1 feature is in good 
agreement with our ATCA observations (see Tabled- The Parkes 
spectrum (Caswell, priv. comm.) shows an absorption feature at 
4—7 km s _1 , which is not present in our data. This absorption 
obscures the red wing of the emission feature and is probably re- 
solved by the ATCA. The Parkes spectrum also contains a hint on 
the —1.1 km s _1 feature at the level of 0.2 Jy. According to our 
measurement (Table|2j the —1.1 km s~ 1 feature is located approxi- 
mately 1 .'6 to the south from the position observed bv lCaswell et alJ 
1 1995). Such an offset implies a beam attenuation factor of 2, as- 
suming the half-power beamwidth of 3'3 at 6.7 GHz. If the flux 
density of the —1.1 km s _1 feature was 0.5 Jy as we had in all ses- 
sions except April 10 (Fig|2j, the attenuate d flux density w ould be 
in agreement with the observations of lCaswell et alj fl995). 

The C425 project from the ATCA archive was devoted to 
search for methanol masers at 6.7 GHz in proplyds in the Orion 
nebula. The observations were carried out on 1995 January 6 and 
March 30 using the 6A and 1.5 A array configurations respec- 
tively. The phase and pointing centre was at 02000 = 5 35 m s , 
#2000 = — 5°23'0", which is about 3'7 away from the positions of 
the 6.7-GHz emission we have measured (Table |2j. The offset of 
3'.7 implies a factor of 2 attenuation due to the primary beam. We 
processed these data the same way as our 6.7-GHz observations 
and created spectra similar to Fig.Q The amplitudes were scaled to 
correct for the primary beam attenuation. Neither the —1.1 km s~ 
feature, nor the 7.2 km s _1 were detected in the January 6 data 
which had a typical la single baseline amplitude rms of about 
0.2 Jy. From the calibrator data we suspect serious technical prob- 
lems during this session. The March 30 data suffer from poor phase 
stability, however, both spectral features are clearly present in the 
data. The 7.2 km s _1 feature has an amplitude of about 0.6±0.2 Jy 
on the shortest baselines, which is in agreement with our measure- 
ment (Table|2| FigQand|2j. The —1,1 km s _1 feature has an am- 
plitude of about 0.2±0. 1 Jy and is revealed only if all baselines are 
averaged together. This value is much lower than the flux density of 
about 0.5 Jy expected from our observations, although it is underes- 
timated because baselines of different lengths have been averaged 
together and the phase stability was bad. However, we cannot ex- 
clude the possibility that the —1.1 km s _1 feature was really fainter 
in 1995, than in 2003, as we have observed it to be varible on a 
timescale of months. 

4.2 The 6.7-GHz emission at 7 8 km s -1 

Inspection of Tables [5] and [5] as well as Figure [3] show that the 
6.7-GHz emission at 7—8 km s _1 is coincident with the 25-GHz 
emission detected in our observations within the uncertainty of the 
measured positions. The tendency to move to the east for lower ve- 
locity exists for the transitions at both frequencies. In Fig. [3] we 



This is a preprint of an article accepted for publication in MNRAS © 2005 RAS 



8 M. A. Voronkov et al. 



overlayed the 25-GHz spots found by Johns ton et aljil993) . which 
have more accurate positions than our measurement but were ob- 
served 15 y ears earlier. Also su perimposed are the 95-GHz image 
obtained bv lWright et alJ ^996) and 3a error ellipses showing the 
positions of the 6.7-GHz emission at 7.2 km s -1 , and 7.5 km s _1 . 
It is clear from Fig.[3]that both components of the 7—8 km s _1 fea- 
ture may be associated with 25-GHz spots, although the agreement 
is better for the 7.2-km s~ component (left ellipse in Fig0 in 
terms of both position and velocity. This component is likely to be 
associated with the south-eastern tip of the crescent shaped region 
traced by the 25-GHz maser spots. The uncertainty ellipse for the 
7.5-km s _1 component may include the southern peak at 95-GHz. 
Althou gh the coarse spectral resolution in the data of )Wright et alJ 
1 1996) makes this comparison difficult. The 7—8 km s _1 feature 
consists of at least two components emanating from close locations. 
In this case the Gaussian fit may suggest the positions of these com- 
ponents to be closer together than they are in reality. Therefore, the 
actual agreement may even be better. The source is too complex to 
determine any definite association with the spatial resolution of our 
6.7-GHz data, but the general coincidence with the area containing 
25-GHz and 95-GHz emission is obvious. A high resolution study 
is required to measure the position of the 6.7-GHz emission with 
respect to individual 25-GHz maser spots. However, it is difficult 
to do this with existing interferometers because the source is very 
weak at 6.7 GHz and most probably resolved. 

It is clear from Fig.0that t he coinc i ndenc e between 25-GHz 
maser spots revealed by Johnston et al. 1 1992) and the peaks of 
the 95-GHz emission detected bv lWright et alJ 17996) is good, but 
not perfect. It seems unlikely that this small discrepancy is due 
to positional errors. Most probably, the 25-GHz masers require 
more specific physical conditions than the 95-GHz maser s. Th e 
bright masers at 25-GHz are relatively rare iMenten et al. 1986), 
while 95-GHz masers are common (e.g.. lVal'tts etal]|200Cl) . An- 
other possible explanation is temporal variability of the 25-GHz 
emission suspected in section l3~2l It is worth noting, that the nature 
of the 95-GHz emission in OMC-1 is not completely clear. This is 
a well known class I maser transition (e.g., Val'tts et al. 2000J) and 
Plamb eck & Wrightl ll988h concluded on the basis of the bright- 
ness temperature that at least a part of the 95-GHz emission in 
OMC- 1 is due to a maser action. How ever, the limit on the bright- 
ness temperature (T b > 500 K) set by Plambeck & Wright! 1 1988) 
is not strong and higher spatial resolution observations are required 
to answer the question unambiguously. A comparison of interfer- 
ometric images of the 25-GHz and 95-GHz masers with that of 
various outflow and dense molecular gas tracers shown that these 
masers arise in the interface regions where outflows interact with 
the a mbient material Ijohnston et alJ 1 19921 : IPlambeck & Mentenl 
1990). Therefore, the 7—8 km s _1 feature is likely to emanate from 
the same volume of gas disturbed by the outflow, even if there is no 
co-existense with the 25-GHz maser spots at finer scales. 

At present we cannot definitively discriminate whether the 
observed 6.7-GHz line peaking at 7.2 km s _1 is a maser or has 
a quasi-thermal nature (i.e., the level populations of the 6.7-GHz 
transition are without inversion). The straightforward approach to 
estimate the brightness temperature limit using the fit result at ei- 
ther velocity listed in Table [2] gives a brightness temperature of 
about 40 K. Because the feature definitely is resolved by our ob- 
servations, this value is a rough estimate of the average brightness 
temperature and should always be considered as a lower limit to the 
actual brightness temperature. 



4.3 Co-existence of masers belonging to different classes 

In the previous section we have shown that the 6.7-GHz (class II 
maser transition) emission feature peaking at 7.2 km s _1 is likely 
to arise from the region where the 25 GHz masers are (class I). The 
spatial resolution of our data is too coarse to be able to claim that 
exactly this is realized in such a complex source as OMC-1, al- 
though there is agreement within the positional uncertainty of our 
measurement. This is a new finding because, as mentioned above, 
methanol masers of different classes co-existing in the same star- 
forming region are usually situated at considerably different loca- 
tions in the source. 

In this section we address the possibility of coincidence of 
the different class maser locations in the strict theoretical sense: 
whether there is a pumping regime allowing a simultaneous inver- 
sion of transitions belonging to different classes in the same volume 
of gas? Previous studies show that this is possible in principle. For 
example the class I II2 — Hi E transition at 26.3 GHz is p resent 
in the l ist of class II methanol maser candidates of Sobol ev"" ;t alJ 
1 1997b). However, the existence of inversion in class II maser tran- 
sitions under conditions where class I masers are strong has not 
been seriously considered, thou gh it was b riefly mentioned in the 
conference proceedings bv lVoronkov et alJ (2005). 

Pumping analysis has shown that the dust temperature of the 
maser environment is one of the key parameters which determines 
the relative strength of the different types of methanol maser transi- 
tions. When the dust temperature is high (T > 120 K) the radiative 
excitation from the ground to the second torsional state becomes 
influential and controls the populations of the masing transitions. 
The levels of the first and the second torsionally excited states are 
located 350 — 650 K and 600—800 K above the levels of th e ground 
state, respectively (Mekhtiev. Godfrey & HougeTJ Il999l . There- 
fore, the corresponding transitions occur at (0.6 — 1.3) x 10 13 Hz 
(20 - 40 A*m) and (1.2 - 1.7) x 10 13 Hz (15 - 25 ^m). The 
strong class II masers appear when the dust temperature is suffi- 
ciently high (T > 150 K) and the gas te mperature is lower o r 
about the same I Sobo lev & Deguchilll994t ISobolev et al]|l997i) . 
Calculations have shown that the strong class I masers do not ap- 
pear when the distribution of population numbers is controlled by 
radiative processes in transitions between the ground and the sec- 
ond torsionally excited state. The situation changes only for rather 
high densities (exceeding 10 8 cm -3 ) which are not likely to occur 
in methanol maser sites. 

The situation at lower temperatures was considered by 
ICragg et'al] ll992l) . This work was limited to the case where con- 
ditions in the masing region are such that the radiative transitions 
between different torsionally excited stat es are not infl uential. This 
implies dust temperatures below 50 K. ICragg et'al] 1 19921) have 
shown that the class I maser transitions are considerably inverted 
only when the gas temperature exceeds that of the ambient radia- 
tion, i.e., excitation processes are dominated by collisions with the 
gas. The class II masers appear only when the radiative excitation 
prevails. So situations where the 6.7-GHz transition can attain con- 
siderable inversion in the star forming region simultaneously with 
the strong class I masers, are limited to cases where the dust tem- 
perature is in the range between 50 K to 120 K and the gas temper- 
ature has a similar value. 

We considered a set of models where the dust responsible for 
the pumping is intermixed with the gas within the molecular cloud. 
The dust temperature was varied over the range 50 K to 260 K. 
Figure 5 displays the results of calculations for models with differ- 
ent dust (Tdust) and gas (T gas ) temperatures, which were assumed 
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Figure 5. The model dependence of the brightness temperature and opti- 
cal depth of various methanol maser transitions on the temperature (T = 
Tgas = Tdust)- Bold lines represent class II transitions. 

to be equal. This is a good approximation for regions with den- 
sity and tempera ture in the range typical for t he 25-GHz masers 
(see, e.g.. ISoboiev. Watson & Okorokovl l200 3 j). which appear af- 
ter th e passage of a moderate shock (e.g., Whitworth & Clarke 
Il997t) . Outlines of the code are des cribed inlSutton et alJ 12004) 
and fu rther detail s can be found in ISobolev et alJ ll997albl) and 
Koooen & Kegel. Il980t). The main difference from the modelling 
procedure described in ISobolevetal1 iT997b) is that the new cal- 
culations take into account radiative processes on the dust which is 
intermixed with the gas within the source, and the external radia- 
tion is reduced to the microwave background with Tt g = 2.7 K. 
The cloud models presented have the following parameters: the hy- 
drogen number density is 10 5 cm -3 , the specific (per unit of the 
linewidth) methanol column density is 10 12 cm~ 3 s, and the beam- 
ing factor (e _1 ) is 10. For all models shown here the optical depth 
of the dust in the source is set to be proportional to the square of 
the frequency and assumed to be equal to 0.5 at 10 13 Hz (30/im). 
The dependence of the brightness temperature and optical depth in 
different transitions versus the temperature (T = T gas — Tdust) is 
shown in Fig[3] A maser line is expected to have a high brightness 
temperature and a negative optical depth. Increasing the tempera- 
ture (T) in such a cloud changes the balance of excitation processes 
in favour of radiative ones. So, the expectation is that the class I 
masers should appear at low temperatures while at high tempera- 
tures the class II masers should shine brightly. The transitions all 
show a single smooth change from non-masing to masing condi- 
tion (positive to negative optical depth), with the notable exception 
of the 6.7-GHz transition. The secondary masing regime of the 6.7- 
GHz transition is at low temperature (60 K), overlapping with the 
regime favouring the class I 25-GHz masers. This is a new pump- 
ing regime, which is not expected in the traditional methanol maser 
classification scheme. 

To investigate this regime further, we repeated calculations for 
a single temperature T = Tdust = T gas = 60 K using a modified 
energy level diagram. We deliberately excluded the energy levels 



Table 4. Optical depths of selected transitions for various models. For all 
models T=60 K. The first 3 models involve the dust intermixed with the 
gas in the cloud (I), the last 3 models involve an external isotropic layer of 
dust (E). 
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Vt = 1 
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yes 
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no 
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6.7 


-15.7 


-16.6 


142 


18.1 


-11.5 


139 


12 


68.3 


68.9 


192 


100 


116 


195 


25 


-12.7 


-12.9 


-10.8 


-11.6 


-11.8 


-10.8 


36 


-11.3 


-11.3 


-10.5 


-10.4 


-10.7 


-10.5 


44 


77.1 


89.6 


-8.9 


21.9 


49.9 


-8.9 


95 


101 


113 


-1.0 


49.8 


72.0 


-1.0 


107 


48.0 


33.6 


192 


49.8 


61.5 


190 



of the torsionally excited states of the methanol molecule to check 
their influence on the overall pumping. In addition, in some mod- 
els we simulated the pumping by an external dust layer as it was 
done bv lSobolev et alJ ll997bl) . For these models an isotropic dis- 
tribution of the dust emission was assumed (i.e. the dilution factor 
is 1). All other parameters were the same as in the models used 
to produce Fig. [5] More detailed analysis require a comprehensive 
investigation of the vast parameter space and will be reported else- 
where. The optical depths of selected transitions obtained in these 
models are listed in Table [4] The first model corresponds to the 
T = 60 K point in Fig.HJ Tabled] shows that the 6.7-GHz transi- 
tion is inverted only if the first torsionally excited state is included. 
Addition of the second torsionally excited state does not siginifi- 
cantly affect the result, except for model 4 involving an external 
dust layer. This means that the pumping of the 6.7-GHz transition 
in this low temperature regime is going through the levels of the 
first torsionally excited state in contrast to the pumping of the bright 
class II masers, which requires both the first and second torsionally 
excited states (Sobolev & Deguchi 1994). In the models where the 
dust is intermixed with gas (labelled I in Table |4j, the 6.7-GHz 
transition is inverted stronger (greater optical depth by absolute 
value) than in the case of the external dust (labelled E in Table |4j- 
In model 4, which includes both the first and the second torsion- 
ally excited states and simulates the pumping by an external dust 
layer, the 6.7-GHz transition is not inverted. Therefore, it is essen- 
tial that the dust is intermixed with the gas for this low temperature 
mechanism to be realized in practice. There is no significant differ- 
ence between the optical depths in the class I maser transitions at 
25 GHz and 36 GHz in these models. This means that these transi- 
tions are pumped mostly through the levels of the ground torsional 
state of the molecule. For the given values of the fixed model pa- 
rameters, the transitions at 44 GHz and 95 GHz are not inverted if 
at least one torsionally excited state is included. These transitions 
belong to the (J + 1) — Ji A + transition series (J=6 and 7, respec- 
tively), which occurs between the same K-ladders as the 5j — 6o A + 
transition at 6.7 GHz, but in the opposite directions. Therefore, the 
tendency for 6.7 GHz to invert due to the pumping through the lev- 
els of the first torsionally excited state reduces the inversion of the 
(J + l)o — Ji A + transitions. Finding the right balance between 
these processes requires a detailed study of the parameter space. 
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4.4 The -1.1 km s" 1 feature at 6.7 GHz 

According to position listed in the Table [2] the —1.1 km s _1 
feature at 6.7 GHz originates in the Orion-South region, 
which contains a prominent dust and NH3 condensation 
(Keen e. Hildebrand & Whitcomblll982l : iBatrla et alJll983l) . How- 
ever, the feature is noticeably blue-shifted with respect to the nor- 
mal spread of velocities for lines tracing this molecular core and no 
emission in methanol lines ha s been previously reported in O rion- 
South at this velocity (see, e.g.. lTatematsu et alll9 93 ; Ment enet alJ 
1988b). The source appears to be unresolved (Tabled which im- 
plies an angular size less than 5" x 2" and a brightness temperature 
greater than about 1400 K. The non-imaging analysis of the visibil- 
ity data described in 13. II yields the size less than 2'.'6 (a fit for all 
data except that from the EW367 array configuration), and thus the 
brightness temperature greater than about 2000 K. This high bright- 
ness temperature and peculiar velocity of the feature indicate that 
the detected line most probably is a maser. The feature appears to 
show temporal variability (see section l3~TV which is another indi- 
cator that it is likely a maser and not quasi-thermal. 

Sensitive observations in the infr ared and radiocontinuum 
IZapata et all2004tlO'Dell & Doil2003l) do not show the presence 
of a compact source at the peak position of the newly detected 
methanol maser listed in the Table |2| However, positional uncer- 
tainties of our observations do not allow us to rule out the possi- 
bility of ass ociation with 1.3 cm sources 139-409 and 140-410 of 
IZapata et alJl2004l) . The nature of these sources is unclear and they 
may be sources of g yrosynchrotron em ission or ionized outflows 
from low-mass stars IZapata et all20o4) . 

It is likely that the methanol maser in Orion-South is associ- 
ated with the outflow. ISchmid-Burgk et alJ il99Ct> detected a CO 
outflow presumably ejected from the CS3/FIR4 infrared source. 
The CS3 and FIR4 positions may actually represent the same ob- 
ject observed in the CS line and 1.3 mm continuum, respectively. 
The nominal separati on between them is 5'.' 5 along the outflow 
axis fSchmid-B urek etal] [l990). The receding part of the out- 
flow is a remarkable 120" long jet with a very narrow width 
(8"). The approaching part of the outflow appears on the other 
side of CS3/FIR4 and has an irregular shape covering the whole 
Orion-South region. The component is seen at velocities from 
to 3 km s -1 and encom passes the methanol maser location. 
Schmid-B urek et alJ <1990h argued that this component of the out- 
flow might protrude into the HII region M42, which could ex- 
plain its ill-defined appearance. Observations of SiO lines trac- 
ing dense parts of the outflow displa y blue wings spa nning down 
to —5 km s _1 IZiurvs. Wilson & MauersbergeJ ll990). The posi- 
tion of the methanol maser at —1.1 km s _1 coincides with the 
blueshifted SiO emission with in the uncertainty of the measure- 
ment (Tableland Ziurys et al.^990). It is therefore possible, that 
the maser is associated with the approaching part of the outflow 
from CS3/FIR4. 

It should be mentioned that the Orion-South region contains 
several outflows (see, e.g.. lO'Dell & D oi 2003) possibly originat- 
ing from the young stellar obje cts whose formatio n was induced 
by outflow descr ibed inlSchm id-Burg k et alJ ll990l). Observations 
of H2 O masers <Gaume et alll998l) and HH objects |5'Dell & Doi 
2003) show that outflows in the region are quite rapid and the ve- 
locity of the newly detected 6.7-GHz methanol maser line compo- 
nent is well within their velocity spread. At present it is not possi- 
ble to associate this methanol maser with any particular outflow or 
young stellar object. Another weak 6.7-GHz methanol maser prob- 
ably associated with an outflow has been found towards NGC 2024 



(source FIR4) bv lMinier etafl [2003). Though NGC 2024 is sit- 
uated rather far from Orion-South (about 4° to the north) it also 

belongs to the Orion complex. 

The excitation analysis of ISobolev et al] ll997al) shows that 
class II methanol masers can be situated at considerable distances 
from the young stellar object if the associated region is affected 
by an outflow which releases methanol from the icy dust grain 
mantles and warms up the dust in the maser vicinity. Great en- 
hancement of the methanol abundance at t he edges of the out- 
flows is confirmed observa tionally (see, e.g.. lBachiller et alll995l 
1 998 . iGibb&Davislf 19981 Association of a considerable portion 
of class II met hanol masers with outflow s is supported by statis- 
tics (see, e.g.. lMalvshev & Sobole-J2003l:ICodella et all2004h and 
observations of pa rticular sources (see, e.g.. lDe Buizer et all2002t 
Minie r et alJ 120031) . We are currently unable to reveal whether 
these weak masers found in the Orion complex are pumped via 
a low temperature mech anism similar to that desc r ibed in 14.31 or 
by the standard scheme iSobolev & Deguchilll994ISobolev et alJ 
1997a b) using the hot dust emission from a nearby source. Ac- 
cording to the position listed in Table|5| the —1.1 km s _1 maser is 
associated with the southern ti p of th e dust condensation in Orion- 
South (Mezger, Wink & Zvlka ll990h . The bulk of the dust in this 
region is rathe r cold (T< 70 K) to allow the pumping via t he stan- 
dard scheme iMezger et alJll990t iMookeriea et al J 12000). How- 
ever, condensations with significantly enhanced dust temperatures 
(Ri200 K) do exist (e.g. JMezger et alll99o\ Due to the lack of the 
dust temperature maps at high spatial resolution, it is impossible to 
discriminate between the pumping regimes on the basis of the dust 
temperature. We detected no 25-GHz emission towards the posi- 
tion of the 6.7-GHz maser in Orion-South. Therefore, if the low 
temperature mechanism is realized, the parameters are unlikely the 
same as used in !4.3l Sensitive observations of the other methanol 
and hydroxyl maser transitions toward 6.7-GHz maser positions in 
Orion-South and NGC 2024 may elucidate physical status of these 
sources. 



5 CONCLUSIONS 

(i) Two 6.7-GHz methanol emission features have been detected 
towards OMC-1, one narrow feature with a peak at —1.1 km s~ 
and a broad feature with a peak at 7.2 km s~ . 

(ii) The 7.2 km s _1 feature consists of at least 2 components 
which are close in both position and velocity. These components 
are likely to be associated with the outflow traced by the 25-GHz 
masers and 95-GHz emission. 

(iii) There is a pumping regime which produces a weak 6.7-GHz 
maser (class II) simultaneously with the 25-GHz masers (class I) 
for relatively low dust temperatures (T ~ 60 K). The pumping oc- 
curs through the levels of the first torsionally excited state and is 
due to the emission of dust intermixed with gas. Brightness tem- 
peratures of the order of 10 6 K can be achieved for such masers. 

(iv) The —1.1 km s _1 feature is most probably a maser with a 
brightness temperature > 1400 K. It originates in the Orion-South 
region and may be associated with the approaching part of the out- 
flow seen in CO. This feature is most probably variable. 

(v) No 52 — 5i E 25-GHz emission has been detected towards 
the Orion-South region at velocities near —1 km s _1 . 
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